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Abstract 

Shock wave formation and propagation in two-dimensional granular materials under vertical 
vibration are studied by digital high speed photography. The steepen density and temperature 
wave fronts form near the plate as granular layer collides with vibrating plate and propagate 
upward through the layer. The temperature front is always in the transition region between the 
upward and downward granular flows. The effects of driving parameters and particle number on 
the shock are also explored. 



PACS numbers: 45.70.Mg, 46.40.-f 



Granular materials are ubiquitous in nature and play an important role in many of our 
industries and daily lives Due to their noncohesive, strongly dissipative properties, 
granular materials behave differently from usual solids, liquids and gases. Under vertical 
vibrations, as the input energy increases, the state of granular materials changes from solid- 
like, to liquid-like and gas-like state. The three states and phase transitions between them 
have been studied experimentally and theoretically m. Fluidized granular materials show 
many interesting phenomena, such as surface pattern Isf, oscillonyi, convection 15|, size and 
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density segregation p, cluster [7|, heap formation j8|| and transport [9|]. Recently sound 
wave propagation attracts much interest because it coexist with most of the above phenom- 
ena and is not well understood. It is found that sound in close packed granular materials 
mainly propagates through force-chains. The speed of sound in sand is approximately the 
same as that in air and very sensitive to the arrangement of grains Computer sim- 

ulations indicate that there exist density and pressure waves in granular materials under 
vertical vibrations Recently time-independent shocks in granular flow past an obstacle 
are observed by experiments and compared with simulation results Moreover the for- 
mation and propagation of shocks in vibrofiuidized granular materials have been studied by 
molecular dynamic simulations and numerical analysis of continuum equations [13]. How- 
ever experimental study of wave propagations in vibrofiuidized granular materials is scarce. 
In this paper we use high speed photography to explore shock formation and propagation 
in vertical vibrated bidimensional granular materials. 

The experiment is conducted with a rectangular container mounted on the vibrating 
exciter (Briiel & Kjaer 4805), which is controlled by a function generator (type HP 3314A). 
We use steel spheres with diameter d = 4 mm and density p = 7900 kg/m^. The container is 
made up of two parallel, / = 90mm(length) and h = 280mm(height), glass plates separated 
by w = 4.1 mm vertically adhered in a Plexiglas bracket. The total particle number N, the 
driving frequency / and the nondimensional acceleration F = An'^f'^A/g {A is the driving 
amplitude and g the gravitational acceleration) are used as control parameters. A high 
speed camera (Redlake MASD MotionScope PCI 2000sc) is used to record the movements 
of the spheres. A frequency multiplication and phase lock circuit is used to generate external 
trigger signals for the camera. The acquisition rate is NpXf with multiple number Np = 25. 
Every recorded image is 22(i(length) x30(i(height). Image processing technique jl4] is used 
to track the locations of all particles. We use vibrating container as reference frame in image 
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FIG. 1: (a) Time-space profiles of density wave Dr{j,k) from 1st to 20th strips, solid circles peak 
points of density at each strip, dashed line polynomial fitting of peak values with time; time = 
corresponds to the time when the plate reaches its maximum height; (b) background density Df, 
versus height, which is below 1 at height above strip 20 (dotted line). Parameters are / = 15Hz, 
r = 5 and = 150. The height of peak point at each strip is set at 90 percent of total height of 
that strip. 

processing. Every image of grains is divided vertically into a number of strips, each with 
width of Id and indicated by 1,2, ... from bottom to top. Assembly average is performed 
over 600 cycles for granular properties at Np phase points in each cycle. Granular density 
D{jd, kAt) is the number of particles in jth strip at time kAt, in which At = l/{Npf). For 
brevity, we omit the units of space and time and use D{j, k) instead of D{jd, kAt) hereafter. 
Granular temperature is defined by T(j, k) = J2n=i — ^b{j, k)\'^/Ns, in which A''^ is the 
number of particles in jth strip at time k, v„ is the velocity vector of the nth particle, and 
the background velocity Vb(j, k) = J2n=i ^n/Ng is the mean particle velocity in jth strip at 
the time k averaged over all cycles recorded. 

As F increases to and beyond 1, granular layer fiuidizes from upper to lower parts and the 
density fluctuations or density wave defined by -Dr(j, k) = D{j, k) — Db{j) (the background 
number density D^^j) = Yl!kLiD{j.,k) /Np) appears in the layer (as shown in Fig. d^a)). 
The bottom of the layer reaches its maximum density at time = — vr/S. Fig. ^a) shows 
the time-space profile of density wave. The fitting curve of density peaks of all the strips 
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indicates that the density wave propagates upward with a nonuniform velocity, faster below 
than above Ath strip and is in agreement with the results obtained from molecular dynamic 



simulations jll|. The wave distorts as it moves upward. It is well known that in ordinary 
gas the distorted wave will evolve into shock wave as it propagates. The calculation of 
Mach number Ma indicates that there exists shock wave in our experiment. As a supersonic 
granular flow encounters an impenetrable plate, velocities of particles (relative to the plate) 
near the plate abruptly decrease to about zero, but those of undisturbed particles are still 
unchanged. This results in a normal shock formation near the plate. Then this shock wave 
leaves plate and propagates upward. We only plot density profiles below 20th strip because 
above that height the number density Dh is too small for accurate statistics (Fig. ^b)). 

In the experiment the shock is identified in a region where the Mach number increases 
from its minimum (< 1) in the disturbed region to the maximum (> 1) in the undisturbed 
region. The Mach number is defined by 



Ma=|^^^^^^|, (1) 

com pom 

HQ, 



c 

where Vbn is the vertical com pon ent of Vb, Vp is the velocity of the plate and c is the speed 
of sound calculated with 



TxdH-xH-^l). (2) 

in which x = 1 + 2(1 + e)z/[l — {i' / I'max)^'^""''^^^]'^ , e is the restitution coefficient, z/(j, k) = 
D{j, k)TTd'^/ (Qlw) is the volume fraction and I'max = Nitd^ / {12hcmiw) = 0.57 is the maximum 
volume fraction JJJ. The height hem of the center of mass at rest is calculated by 

Km = ^[{l - V3/2)n, + {V3/2)nl] + ^{1 + VSn,), (3) 
in which = l/d — 0.5 is the average number of particles per strip, = int(N/nh) is the 



number of full strip and tiq = N — nhUh is the particle number at the highest strip [1^. We 
investigate the propagation of the shock by examining the distribution of granular density, 
temperature and Mach number at four time in one cycle (shown in Fig. |2I). The maximum 
Mach number in t 
C = lw/{V87idD) 



le shock and the shock width relative to the mean free path of particles 
l^l as a function of time are drawn in Fig. 121 to show the shock dynamics. 
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FIG. 2: Velocity field (left column), Mach number Ma(asterisk) and scaled temperature(open 
circle) (right column) as a function of lieight(right) at four different times. At time = vr the plate 
is at its lowest position. Parameters here are the same as those in Fig^ 
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At time = 1.20n, the plate moves upward with Vp = 0.306m/s and begins to colhde with 
the granular layer. Particles at lower strips begin to be compressed by collision. In the 
mean time most of the particles in the undisturbed region are still falling towards the plate 
with supersonic speed. This results in the formation of a shock wave. In the shock region 
the Mach number increases with height and reaches a maximum value 2.73 at the 5th strip. 
The particle velocities are randomized by collision and the granular temperature decreases 
with height in this region. 

At time = l.GOvr, the plate moves upward with Vp = 0.495m/s. The granular layer 
continues to be compressed on the plate, more and more particles come into the compressed 
region, resulting in the propagation of the shock up through the layer. In the shock region, 
the Mach number increases with height and reaches its maximum (Ma = 2.95) at 10th strip 
(Fig. 3). A temperature peak appears at 3rd strip which is the bottom of shock region and 
is the transition region between upward and downward granular flows. 

At time = 2.007r, the plate reaches its maximum height. The granular layer leaves the 
plate and flies freely. The shock propagates upward with the maximum Mach number 
decreasing to 1.39 at the 12th strip. As seen in Fig. 13 the shock width now decreases to 
its minimum (approximate 1.7 times the mean free path), indicating that the shock is fully 
developed. The steepen temperature wave front propagates upward to 10th strip. 

After time = 2.007r, the plate begins to move downward. The maximum Mach number 
in the shock region decreases. At time = 2.167r, the shock disappears because the maximum 
Mach number at strip 13 damps to be less than unity. The temperature wave front continues 
to propagate upward until disappears in the dilute region of the layer. Then the granular 
layer is ready for the next collision with the plate. 

To investigate how changing F and N affect the shock, we perform the experiment with 
F from 5 to 8 and from 60 to 250. As F increases, the velocity of granular layer relative 
to the plate increases, resulting in the increase of maximum Mach number and the speed of 
shock. On the other hand, as F increases, the average number density decreases (particles 
expand to more space of the container), and this results in a decrease of the propagating 
velocity of the wave [l^. As a result the total effects of these two opposite factors make a 
little change of the velocity of shock with the change of F (Fig. EJa)). As total number of 
particles increases from 150 to 250 (the layer depth from 7d to lid) with driving parameters 
fixed, the average number density increases. As indicated in Fig. |^b), denser layer collides 
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FIG. 3: The peak value of Mach number and shock width Wg/C relative to the mean free path 
as a function of time. The shock begins to form at time= 1.207r and disappear at 2.007r. The 
shock width decreases to about 1.5 at time= 2.007r, indicating that the shock is fully developed. 
Parameters here are the same as those in Fig. ^ 

with the plate later than more dilute layer. Thus increasing layer depth causes the shock to 
form later in the cycle and the shock propagates faster through the higher density layer. 

In conclusion, we prove the existence of shock wave in vibrofuidized granular materials 
experimentally. The shock forms as particles collide with vibrating plate and propagates 
upward with a steepen temperature front in the transition region between upward and 
downward granular flows. The velocity of shock depends on the velocity of the plate when it 
collides with granular layer and the number density of granular layer. This is helpful for the 
understanding of phenomena in vibrofluidized granular materials such as surface instability, 
convection and energy transfer, etc. 
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